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ABSTRACT 

Recent WMAP observations indicate that the primordial abundance of hthium should 
have been 2-4 times higher than the value measured in low-metallicity stars. We show 
that this discrepancy can be explained by the process of ambipolar diffusion in pro- 
tostellar clouds, which leads to depletion of elements with low ionization potential, 
such as lithium and sodium. In high-metallicity stars, [Fe/H] ^ — 1, the depletion of 
lithium due to this process is found to be negligible, but for [Fe/H] ;$ — 1.5 the lithiimr 
abundance falls by a factor of ^ 3, consistent with WMAP observations. 
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1 INTRODUCTION 

Since the discovery of a nearly constant lithium abundance 
(Spite & Spite 1982) in the metal poor stars, there has been 
a tendency to assume that the measured value reflects the 
primordial abundance. However, in recent years this inter- 
pretation came into conflict with the high value of f2b ob- 
tained by the WMAP satellite (Spergel et al. 2003; Spergel 
et al. 2006), which in the standard big bang nucleosynthesis 
model corresponds to a factor of 2-4 larger lithium abun- 
dance. In this paper we show that the process of ambipolar 
diffusion in protostellar clouds can provide a natural expla- 
nation to the missing lithium problem. 

The process of ambipolar diffusion, i.e., an outflow of 
charged particles from neutral gas caused by magnetic force, 
is primarily thought of as a way to reduce the strength of 
magnetic field. It has been long recognized that in addi- 
tion to efficient cooling, which removes the thermal pressure 
support, the collapse of the diffuse gas into a star requires 
some mechanism to reduce the magnetic pressure, for which 
purpose ambipolar diffusion seems to be the most promis- 
ing candidate (e.g. Mestel & Spitzer 1956; Nakano 1976; 
Mouschovias & Paleologou 1981; Mouschovias et al. 2006). 
As we show in this paper, besides reducing the magnetic flux, 
ambipolar diffusion can also affect the chemical composition 
of the gas by preferentially expelling elements with low ion- 
ization potential. In particular, we flnd that in the metal- 
poor stars ([Fe/H] ^ — 1.5) * the abundance of lithium falls 
by a factor of ~ 3, which allows to explain the discrepancy 
with the WMAP resuhs. 

The paper is organized as follows. In §2, we discuss how 
ionization balance evolves in a collapsing protostellar cloud. 



[X/H]= Logw{n{X)/n(H)) - Logw(n(X) /n{H))Q 



In §3, we estimate how the chemical composition of the cloud 
changes during the collapse. In §4, we summarize our results. 



2 IONIZATION BALANCE 

Partial ionization of protostellar clouds can be sustained by 
several different processes. Ionizations by cosmic rays, which 
are believed to be dominant at the present epoch, initially 
produce mostly hydrogen and helium ions. However, since 
the charge-transfer reactions are typically much faster than 
radiative recombinations, the light elements pass the charge 
to metals with lower ionization potential, which therefore 
make up most of the ions (Oppenheimer & Dalgarno 1974). 
In equilibrium the ionizations by cosmic-rays are balanced 
by recombinations 

C,n = ax^ri" , (1) 

where C, is the cosmic rays intensity, a is the recombination 
coefficient, n is the gas particles number density and x is 
the fractional ionization. Neglecting the possible change in 
a, which may result from change in temperature and chem- 
ical composition of the gas, we flnd from eq. (1) that the 
fractional ionization scales with density as a:: oc n~^^'^ . 

At high redshifts, when the UV flux was significantly 
higher and the cosmic-ray flux may have been lower, ion- 
izations by UV photon with energies below 13.6 eV (hy- 
drogen ionization potential) to which the neutral hydrogen 
is transparent, may have been the dominant source of ions 
in protostellar clouds. In the optically thin limit x is again 
proportional to n~^l'^ . At high densities the dependence of x 
on n becomes non-trivial, as metals with high ionization po- 
tential become self-shielded, while those with low ionization 
potential (which can be ionized by wider spectral range) are 
still optically thin. Nevertheless, up to the point when all 
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abundant metals become self-shielded, the dependence of x 
on n should be only slightly steeper than in the optically 
thin limit. 



3 PROTOSTELLAR CLOUD COLLAPSE 

We consider the collapse of a protostellar cloud under its 
own gravity, starting from a point of magnetostatic equi- 
librium. Following previous authors (e.g. Chandrasokliar & 
Fermi 1953), we ignore the detailed structure of a cloud, as- 
suming instead that it can be described as an infinite cylin- 
der with a single set of characteristic values of density, mag- 
netic field strength and ionized fraction ^ . 

The rate of contraction of a gas cloud is determined by 
the competing forces of gravity and magnetic pressure 



dt'^ GTrpr -|- 2^pj. ' 



(2) 



where r and p arc, respectively, the characteristic radius and 
density of the cloud. The magnetic field evolves according 
to 

^ = Vx(«ixB) = -^B, (3) 

where Vi is the velocity of the charged particles, which ex- 
ceeds the velocity of the neutrals by 
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■ Vi — 
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(V X B) X S 



(4) 



where 70 « 9 x 10^^ cm^s^^g"^ (Osterbrock 1961; Draine 
et al. 1983). The abundances of metals evolve according to 

6v 



-( 

dt V 



dt \ riH / 



rim 
Tin 



(5) 



where Xm is the metal ionized fraction, ExmT^m = xn. 

We assume that initially the relative abundances of met- 
als are solar. Further we assume that x scales with density as 
a power law, x oc p'^ , and that all ions are made up by metals 
with lowest ionization potential. Setting the time-scale for 
collapse of protostellar clouds, tcou, which is generally esti- 
mated to be of order 10 Myrs (e.g. Mouschovias et al. 2006), 
allows us to chose the initial value of x, which is almost di- 
rectly proportional to tcoU- Using these initial conditions, we 
integrate the equations (2)-(5) for different mctallicities. 

The results of the integration for lithium, sodium and 
iron abundances are illustrated by Figures 1-3. At high 
mctallicities the impact of ambipolar difi^usion on element 
abundances is negligible, eis the number of ions expelled by 
magnetic field is too small compared to the number of metal 
atoms. Increasing tcoii, which allows to increase the number 
of expelled ions, or decreasing metallicity can make the im- 
pact of ambipolar diffusion more significant, but the deple- 
tion does not generally exceed a factor of a few. This limi- 
tation results from the fact that when magnetic flux drops 
by a factor of a few (with the corresponding drop in metal 
ion abundances) it can no longer effectively slow down the 
gravitational collapse and the dynamical time begins to fall 



t We have verified that our results are almost Insensitive to the 
chosen geometry by repeating the calculations for a spherically 
symmetric cloud. 




Figure 1. Lithium depletion as a function of stellar iron abun- 
dance. The solid and dashed lines arc for x oc ri~"'^ with t^oii = 10 
and 5 Myrs, respectively. The dot-dashed and dotted lines are for 
X oc n"*^'® with t^^on = 10 and 5 Myrs, respectively. 




Figure 2. Same as Fig. 1 for sodium. 
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Figure 3. Same as Fig. 1 for iron. 
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much faster than ambipolar diffusion time-scalo. Thus, at 
low metallicities the impact of ambipolar diffusion is only 
weakly sensitive to the assumed parameters (/3 and tcoii), 
while at high metallicities it is generally negligible. On the 
other hand, the mctallicity at which the transition between 
the low and high metallicity regimes occurs is very sensitive 
to tcoll- 

Several uncertainties in our model should be pointed 
out. So far we have assumed that all metals are in the 
atomic phase, neglecting the depletion on dust, which is 
highly uncertain for low metallicity/high redshift regions 
(e.g. Maiolino et al. 2006). Since most of the dust parti- 
cles are charged, they can be also be expelled from gas by 
magnetic pressure, which opens an additional way of reduc- 
ing metal abundances (Ciolek & Mouschovias 1997). An- 
other uncertainty is the magnitude of the magnetic field at 
high redshifts. Observations indicate existence of a strong 
magnetic field already at z > 2 (Athreya et al. 1998). How- 
ever, at even higher redshifts stars might be forming in gas 
clouds with dynamically weak magnetic fields. Recently, Silk 
& Langer (2006) have shown that even weak magnetic field 
would be sufficiently amplified during cloud collapse to be- 
come dynamically important, but in such a scenario the im- 
pact of ambipolar diffusion on metal abundances is very hard 
to estimate. 
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4 SUMMARY 

We have shown that in metal-poor clouds with dynami- 
cally important magnetic field the process of ambipolar dif- 
fusion results in lithium depletion by a factor of a few, 
which allows to reconcile WMAP measurements with the ob- 
served abundances. Likewise, the depletion of sodium may 
explain why [Na/Fe] ratio is constant in high metallicity 
stars ([Fe/H] ^ — 1), where the impact of ambipolar diffu- 
sion on the abundances is weak, but drops in the metal-poor 
ones (Carretta et al. 2000). It should be noted, however, that 
since, unlike lithium, metals (such as Fc and Na) are not of 
primordial origin, their abundance patterns may have other 
causes besides diffusion. 

In the metal-rich clouds the impact of ambipolar diffu- 
sion on the lithium abundance is negligible. However, by 
reducing the lithium abundance in the metal-poor stars, 
ambipolar diffusion is increasing its abundance in the inter- 
stellar medium. Therefore, later (i.e. metal-rich) generations 
of stars may in fact form with increased lithium abundance. 
Since the magnitude of this increase depends strongly on the 
previous star-formation history and the inter-stellar gas dy- 
namics, there should be a significant scatter in the lithium 
abundance of the metal-rich stars, which is consistent with 
observations. 



ACKNOWLEDGMENT 

This work was supported by the W.J. McDonald Fellowship 
of the McDonald Observatory. 



© 0000 RAS, MNRAS 000, ??-?? 



